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Abstract 

We propose new canonical form of the fluctuational Hamiltonian which takes 
into account the fact that in the vicinity of the critical point there are two 
fluctuating fields. They are the field of the number density and the entropy. 
The proposed canonical form is based on the catastrophe. In contrast to 
the standard approach of Landau- Ginzburg Hamiltonian which is based on 
^3 catastrophe the canonical form proposed for the fluctuational Hamiltonian 
naturally includes the asymmetric coupling between the fields. 



1. Introduction 

The modern theory of the critical phenomena in liquids is based on the 
Landau- Ginsburg fluctuation Hamiltonian: 



HLcVPir)] = Jdr 



(1) 



where all designations are standard and the field ip is normalized so that the 
coefficient at the square gradient is 1 fl], 3 . Such a form of the fluctuational 
Hamiltonian is sufficient to describe the critical asymptotical behavior of the 
free energy and is the basis for the isomorphism of such behavior with that of 
the Ising model. Nevertheless, the Hamiltonian lacks the asymmetrical terms 
which in general appear in the description of the molecular fluids 0, M. 
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The asymmetry of the coexistence curve (CC) is the characteristic feature 
for one-component hquids and solutions. It is absent in the Ising model and 
numerous lattice models. This circumstance indicates that the main role 
in the appearance of the asymmetry belongs to the effects of finite size for 
molecules of liquids. 

The simplest characteristics of the asymmetry degree for the CC is its 
diameter defined as 

1 

'fd = -{^i + Vv), (2) 

where is the order parameter and subscripts /, g stand for the coexisting 
liquid and gas phases correspondingly. The latter can be identified with the 
extensive variable suitable to differentiate between the coexisitng liquid and 
vapor phases. As rule, the main attention in literature is paid to the diameter 
of the CC in terms density-temperature. Its investigation has a long history. 

In Pokrovskiy's work ||^ (see also 0) it was shown that the classical law 
of the rectilinear diameter ^ is violated: 

n, = - 1 = /^i-.kr-" + Di|r| + . . . , (3) 

where n is the number density, a is the critical exponent for the heat capacity. 
The new term in Eq. (j^) is caused by the fluctuation effects. In accordance 
with [H this contribution is connected with the non-orthogonality of the 
fluctuations of the density and the energy. From another point of view it 
means that the expansion of the density fluctuations with respect to the 
algebra of strongly fluctuating variables takes the form (see [0): 

fi = Ai + \A2 + ... (4) 
Since ^ ^ = ±ai|r|^ (different signs correspond to the liquid and vapor 
branches of the CC) and I A2) = a2|T|^^", the expansion Eq. (H) directly to 



Eq. (^. For laboratory order parameter, which in general is nonsymmetrical, 
the experimental results should be fitted by the combination P, M: 



(Pd = D2^\t\'" + Di_„|rr" + Di|r| + . . . , (5) 

where /3 is the critical exponent for the CC in terms temperature-density. The 
appearance of such singularity in the diameter is trivial if is a nonlinear 
function of some initial (laboratory) order parameter. The appearance of 



2 



the r^^-singularity for the number density is partly surprising. The density 
is the statistical average of the simple microscopic field of the microscopic 
density and is nothing but the unary distribution function. This means that 
the density can be considered as the simplest irreducible quantity in case 
of the liquid-vapor CP. So if the 2/3 anomaly exists for this quantity then 
the question about the nature of r^^-singularity of the density is deeper and 
does not reduce to trivial algebraic relations between the thermodynamic 
averages. 

The revival of interest to the problem of |rp^- term is connected with the 
complete scaling proposed in This approach is based on the supposition 
that all three thermodynamic fields temperature r, the chemical potential 
/i = (/i — /ic)//^c and the pressure p = {P — Pc)/Pc are mixed equally. Then 
it is easy to show that n is a nonlinear function of the average of the scaling 
fields which represent two strongly fluctuating quantities. As a result, the 
appearance of the \t\'^^- term in Eq. is caused by the nonlinear algebraic 
dependence between n and the average of the order parameter. Numerous 



application of this approach are given in |T2|-p!5|. In particular, in complete 



scaling approach the order parameter for the molecular liquids is the non- 



linear combination of the density field and the density of the entropy |T6 
The density of the entropy is weekly fluctuating field but the density field 
fluctuates more strongly. 

From another point of view the problem of the singular terms in Eq. 
is solved in |T^, ^ on the basis of so called the canonical scaling. The last 
is the synthesis of the scaling approach and the canonical formalism, devel- 
oped in the catastrophe theory [^. In these works the regular procedure 
for the construction of the fluctuation Hamiltonian for liquids, which has the 
Ising-like form, is developed. Note, that the Ising-like Hamiltonian coincides 
with that given by the standard canonical form (see As will be shown 

below it corresponds to the thermodynamic potential with the only one van- 
ishing eigenvalue of the stability matrix for the Gibbs potential. The CC 
for liquids and solutions in the canonical variables is fully symmetrical. One 
can say that the canonical order parameter can be considered as the optimal 
one. The asymmetry of the CC appears only if one returns to the labora- 
tory variables. The evident relations between the canonical and laboratory 
order parameters as well as between the canonical and laboratory fields are 
established. The distinctive peculiarity of these relations is their nonlinear 
character. Nevertheless, in contrast to the complete scaling approach the 
canonical one operates with two fluctuating fields and uses the nonlinear 
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transformation rather for the field variable of the order parameter than the 
nonlinear transformation for its average value [|^. 

The motivation of this work is the search for the unifying Hamiltonian 
formulation of the problem of the asymmetry suitable for the determination 
of the proper order parameter in case of two fluctuating fields one of which 
fluctuates less than the other. 



2. Gibbs potential analysis 

Now we discuss briefly how the order parameter ip and the field h con- 
jugated to it, as well as the value of r in Eq. (|l]) are connected with the 
expansion of the fiuctuational deviation for the Gibbs potential g{n^ s) near 
the critical point: 

(6) 



Here 



'^2,0 



dn 



dT/dsl 



dT 
ds 



T 



(7) 



!f = n-nc, C = s - Sc + j; dT/dn\ 



in - n. 



where n is the chemical potential, c„ is the heat capacity, n and s are the 
density and entropy correspondingly, the subscript "c" denotes their values 
at the critical point. 

Usually it is taken into account that only one of two coefficients 02,0 and 
ao,2 in (0) vanishes at the critical point. The supposition 



(^2,0 = -TT- 
on 



^ , T ^Tr 



(8) 



corresponds to the situation when only the fluctuations of the density be- 
come anomalously large. Note that in the framework of the thermodynamic 
fluctuation theory the second coefficient 00,2 7^ 0. Indeed, in the vapor phase 
Cy takes the finite value and the transition to liquid state is accompanied by 
the addition of the only degree of freedom (the order parameter), which leads 
to the increase of c„ on ^/c^. 
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Since ao,2 7^ 0, the fluctuations of ( are bounded and we can neglect them 
thus reducing to the standard Landau expansion: 

1111 
SgL = 2«2,oV5^ + 2«3,oV5^ + ^«4,oV5'^ + -^a5,of^ + ... . (9) 

In fact the local part of Eq. (P is nothing but the local form of Eq. (Q) 
truncated to the minimal polynom sufficient to describe the reconstruction 
of the minima which represent the coexisting phases. 

The transformation from the infinite series (|^) to the canonical form (|ID is 



discussed in details in In that case when Cy is bounded, the canonical 



form for the local part of the fiuctuation Hamiltonian takes the form: 

hiocW = SqlI^^^^^) = ^gij^ + ^r^^ + h^ip . (10) 
The canonical order parameter ijj is the nonlinear function of initial variable 

1 2 1 3 

= <^ + 2 ^aV' + ■■■ ■ 

Explicit expressions for g,r and as well as for the coefficients 7„c an be 
found in . 



However, in many cases the applicability region of the thermodynamic 
fluctuation theory or the mean-fleld approximation is absent. Argon can 
serve as the simplest example of such a kind. It is evident from the behavior 
of its isothermal compressibility. In other words the Ginsburg temperature 
for argon and many other liquids is close to unity TGi{Ar) ~ (0.1 -^ 0.3) 
p2| . It means that the behavior of the heat capacity is determined by 
fluctuation effects. Therefore, fluctuation effects cannot be ignored and they 
lead to the vanishing of the quantity T/c„, which is the second eigenvalue 
of the quadratic form in Eq. (U). In such a case the canonical form for the 
local part of the fluctuation Hamiltonian should be generalized since both 
variables ip and ( should be considered as the strongly fluctuating ones. 

The new version of the fluctuational Hamiltonian will be able to reflect 
more exactly the relative role of mechanical and thermal variables, which are 
used for the description of the critical behavior. At that, if an order parame- 
ter is modeled by a mechanical variable, the auxiliary fleld in the fluctuational 
Hamiltonian is given by caloric one and vice versa. In particular, the specific 
volume as the order parameter is completed by the specific entropy. Tak- 
ing into account of this variable is especially important for complex liquids. 
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molecules of which rotate, can form H-bonds and have internal degrees of 
freedom excited in the thermal motion, as well as for multicomponent solu- 
tions. Since fluctuations of the specific volume and specific entropy are not 
orthogonal, the increase of the role of caloric effects should be accompanied 
by a growth of the coexistence curve asymmetry. In connection with the last 
problem, it is necessary to note that the asymmetry of coexistence curves 
arises not only due to caloric effects. Another mechanism is caused by hard 
core effects. Their manifestation in complex liquids is similar to that for 
simple liquids. Therefore in this paper our attention will be mainly focused 
on the inclusion of the caloric effects into the fiuctuational Hamiltonian. 

Therefore the fiuctuational Hamiltonian should be built on the basis of 
the canonical form for the thermodynamic potential which is consistent with 
the fact of the existence of two strongly fluctuating fields. Another support 
for such conjecture is given by the two-dimensional lattice models where 
the requirement of the conformal symmetry leads to the two independent 



strongly fluctuating fields From here and the very general isomor- 

phism principle it follows that the bare fluctuation Hamiltonian for liquids 
should also contain two independent order parameters. 

3. Fluctuation Hamiltonian of the system 

The local part of the fluctuation Hamiltonian for a system, for which two 
eigenvalues 02,0 and ao,2 of the stability matrix tend to zero at the critical 
point, in accordance with the catastrophe theory must be built on the basis 
of catastrophe [p] : 



^9\^^i,{^,0,C^<^(v,0 = +-r^lj +-ua +h^^ + Ka, (11) 

where ip and a are the corresponding canonical variables which are the non- 
linear functions of the initial variables if and (. In view of the results of 



Ref. |18[ the parameter of the asymmetry coupling b canonically represents 
the degree of the asymmetry of the initial Hamiltonian. 
The fiuctuational Hamiltonian consists of two parts: 

= Hloc + Hql 

OC ; 

the local part if/oc 



Hloc = j Sg{i^{r), a(r)) dr 
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and the quasilocal part Hgioc- For the latter we use common quadratic gra- 
dient approximation. In accordance with Eq. ([TT|) the local part of the fluc- 
tuational Hamiltonian has the form: 



Hi,Mr), air)] = Hi^M^)] + H^oMr)] + H^Hr), ^(r)] , (12) 



with 



dr 



HioMr)] -- 
HioMr)] = j dr 
Hint[ip{r),a{r)] = j drbipa 



1 -.2 , 1 .,,4 



-h^ip + -rip + -gi^ 



U I 2 



(13) 
(14) 
(15) 



Here ip{r) and a{r) are the fluctuating fields. In accordance with Eq. (P) tp{r) 
can be identified with the density while a{r) is the combination of the density 
and the entropy fluctuations. The local part of the Hamiltonian should be 
completed with the relevant gradient terms. They are 



^ (V^(r))'rfr , 



Hniocl-ipir)] = 
HnioMr)] = / ^ V^(r) ■ Va(r) dr . 



(16) 
(17) 



We did not include the quadratic term (Va)^ because as we show below 
within the dimensional analysis it is irrelevant. 

3.1. Simple dimensional analysis 

Here we perform the analysis of the Hamiltonian (|1^) introduced above 
using the simple scaling [^. From the elementary analysis of dimensions 
and the fact that the main contribution to the coefficient u are caused by 
the fluctuations of the free field ip it follows that its dimension is Ai°^ = 
A^°^ = A — d = e. The character of the temperature dependence of r does not 
change in comparison with (|l|), so Ar°^ = 2 and r = r^^r. Also it is obvious 
that Act°^ = 2A^'''' = d — 2. Thus for the conjugated field we have A||^'' = 

d — Ai°^ = 2, thus it is nothing but the canonical temperature variable. The 
field obviously plays the role of the canonical chemical potential and the 
condition = determines the coexistence line. The scaling dimension of b 
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- the strength of asymmetric interaction is A[°'' = | ( y ~ ) ? which becomes 
relevant below d = ^ ^ind thus equivalent to the quintic asymmetrical 
term in the standard consideration The gradient term in Eq. ( pTD has 
the same dimension as the asymmetric term (V^/')^, which should be taken 
into account along with the quintic coupling [|2^. The term (Vcr)^ can be 
omitted because it is irrelevant (with dimension Ay^ = 2 — d = —2 + e). 

For the Ising model because of the symmetry — )■ —ip, b = c^cr = and 
the fields ip and a decouple. 

Previous consideration of the classical scaling dimensions obviously allows 
to identify a with the field : ^/'^(x) :. With account of the fluctuations the 
dimensions are renormalized = A^'* + ^ and A^. = A^j^ + 2 — ^ so 
that = \jv and indeed /i^ is a temperature field. Here v is the critical 
exponent of the correlation length and ri is the Fisher's critical exponent of 
the anomalous dimension ||28| . The dimension of the parameter u is A^ — ^ 



V 



d = ^. This leads to the conclusion that indeed u oc r", i.e. asymptotically 
it is the inverse value of the specific heat. 



4. Fluctuation Hamiltonian for binary solution 

In this Section we briefly discuss the specificity of the fluctuation Hamilto- 
nian in binary solutions. The method of consideration is naturally generalized 
for arbitrary mult i- component solutions. 

Let us consider the binary solution, for which the total number of particles 
remains to be constant. In this case the thermodynamic state of the system is 
set by the temperature and pressure as well as the difference of the chemical 
potentials of components: /i = /i2— yUi- Considering yU as a function of specific 
volume V, entropy and concentration x we can write the following expression 
for the increment of the chemical potential near its minimal value: 



^Sx^ + ^{5x + c5vf + Mi{5x + a5s + h5vf 



+ ... 



:i8) 



Here Mj, i = 1,2, 3, are the main minors of the matrix: 



dT/ds dT/dv dT/dx 
M= I dP/ds dP/dv dP/dx 
dfi/ds dfi/dv dfi/dx 
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In fact, the diagonal form ( [T^ ) is naturally obtained by the Lagrange 
method. 

The coefficients a, b and c are some combinations of the elements of M. 
It is not difficult to verify that near the mixing critical point: 



M3 
Mo 



dx 



P,T 



A ^2 T 

and — — ~ — 

Ml c„ 



Since Mi = — and Cxv takes the finite value for binary solution the 
fluctuations of the combination Su = Sx + aSs + bSv are bounded and the 
contribution of corresponding term in Eq. ( [T8| ) can be omitted. 

Let us compare (|l^) with the analogous increment of the chemical po- 
tential for one-component liquid near its vapor-liquid critical point. In this 
case 



A2 1 
fJ, = - 



Ml ^ ^ 



+ 



where 



M 



dT/ds 
dP/ds 



dT/dv 
dP/dv 



8P_\ 

dv It 



and Ml 



As we see, the characters of the near-critical behavior of 5^yU for binary 
solutions and one component liquids are fully identical. Therefore we con- 
clude that critical phenomena in them are isomorphic, i.e. their fluctuation 
Hamiltonians have the same form. 



5. Conclusions 

We have performed the stability analysis of the thermodynamic potential 
for the molecular systems with the critical point where tgi is not small ^ 
0.1. Such situation means that there is no crossover region to the mean-fleld 
behavior so the canonical form for the fluctuational Hamiltonian include 
two independent fluctuating fields. Both fields are strongly fiuctuating but 
one of them has the strongest singularity caused by the slowest decay of its 
correlator. The proposed canonical form of the fiuctuational Hamiltonian is 
constructed on the basis of catastrophe and contain asymmetrical term 
which describes the coupling between these fields. If the contribution of less 
fluctuating fleld is neglected we get the standard Landau-Ginzburg- Wilson 
Hamiltonian. The facts that is not small and the leading correction to 
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the quadratic part (random phase approximation) is determined by ip^ term 
allow to conclude about global cubic character of the binodal. This is in 
accordance with the global isomorphism proposed in Ref. . 

Simple scaling analysis shows that one of the field can be identified with 
the density-like one while another field corresponds to the fluctuations of the 
entropy. We hope that the proposed canonical form can serve as the basic 
fluctuational Hamiltonian which describes the asymmetry effects character- 
istic for the real fluids more adequately. The problem of derivation of (0) on 
the basis of microscopic Hamiltonian can be performed within the canonical 
approach [0, ^ . 

The physical reasonings for the appearance of two fluctuating flelds are 
connected with the account of the hard core effects in the microscopic Hamil- 
tonian. Usually the microscopic Hamiltonian is divided into short-range and 
long-range parts due to the corresponding division of the interaction potential 
into repulsive and attractive ones. The repulsive part is usually integrated 
over and the fluctuational part depending only on the density fluctuations 
appears (see e.g. Refs. [|], 3). The asymmetry of the coexistence curves for 
real liquids arises due to the breaking of the particle-hole symmetry for the 
lattice gas. This symmetry is flrst of all violated because of different role 
of the hard core effects in liquid and vapor phases. It is also violated by 
different manifestation of caloric effects in coexistence phases. 

It is necessary to stress that the canonical fluctuation Hamiltonian (|12D 
of liquids is not isomorphic to that for the Lattice Gas or equivalently for 
the Ising model, since the interaction Hamiltonian (^) is an odd function of 
the order parameter: Hint[—ip{r),cr{r)] = —Hint[ip{r),(7{r)]. We can speak 
about isomorphism between them only if Hint[ip{r),a{r)] is negligibly small. 



The results of Ref. [|I8| allow to connect the parameter of the asymmetry 
coupling b with the asymmetry of the initial Hamiltonian on the basis of 
the model for the equation of state for the real molecular fluids. With such 
model equation of state the caloric and the hard core effects can be taken 
into account explicitly. The detailed consideration of this question will be a 
subject of prospective work. 
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